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Abstract: The Southern part of the Pierre Auger Observatory is nearing completion, and has 
been in stable operation since January 2004 while it has grown in size. The large sample of data 
collected so far has led to a significant improvement in the measurement of the energy spectrum 
of UHE cosmic rays over that previously reported by the Pierre Auger Observatory, both in 
statistics and in systematic uncertainties. We summarize two measurements of the energy 
spectrum, one based on the high-statistics surface detector data, and the other based on the 
hybrid data, where the precision of the fluorescence measurements is enhanced by additional 
information from the surface array. The complementarity of the two approaches is emphasized 
and results are compared. Possible astrophysical implications of our measurements, and in 
particular the presence of spectral features, are discussed. 



UHE cosmic ray energy spectrum 

The Pierre Auger Observatory measures ex- 
tensive air showers induced by the highest en- 
ergy events [E > 10^^ eV) using two detection 
techniques. Firstly, a collection of telescopes 
is used to measure the ultraviolet fluorescence 
light produced when electrons in the shower 
excite nitrogen molecules in the atmosphere. 
This technique will be referred as FD (Fluores- 
cence Detector). It measures the longitudinal 
development of the air-shower and can only be 
used during dark and moonless nights, yielding 
a duty cycle of roughly 10%. The second tech- 
nique (called SD for Surface Detector) uses an 
array of water Cherenkov detectors to sample 
the shower front at ground level. The SD has a 
duty cycle of 100% and the detection efficiency 
is 100% for energies above lO^^-^ eV (lO^^-^ eV) 
at zenith angles below (above) 60°. The show- 
ers recorded by the SD are quantified in size 
using the reconstructed signal at 1000 m from 
the shower axis, called S(IOOO) [4]. At large 
zenith angles (above 60°), due to deflection of 
the shower particles in the geomagnetic field, 
another energy estimator iVig is used [2]. The 



conversion from these two SD estimators to the 
primary energy could be calculated using full 
Monte Carlo simulations but the lack of knowl- 
edge of the primary mass and the uncertainties 
in the hadronic models introduce large system- 
atics. Therefore we use a subset of showers 
called hybrid events that are detected by both 
the SD and the FD. The conversion parameters 
from the SD estimators to the energy measured 
by the FD then are derived experimentally. A 
comparison of the results of this calibration 
with the expectations from Monte Carlo sim- 
ulation can be found in [8]. The FD measures 
fiuorescence light in proportion to the energy 
deposited by the shower, and so the technique 
is calorimetric. There is, however, a small cor- 
rection to account for the energy deposited in 
the ground by high energy muons and neutri- 
nos. This "invisible energy" correction has a 
small dependence on mass and hadronic model. 
The applied correction is based on the average 
for proton and iron showers from the QGSJet 
model. This correction factor is about 10% 
and its systematic uncertainty contributes 4% 
to the total uncertainty in FD energy [17, 10]. 
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the relative normalization factors necessary to 
match the spectra with each other. We find 
that the different spectra are in excellent agree- 
ment with normalization factors smaller than 
3%. We combine the three spectra weighting 
each bin based on its statistical uncertainty. 
The final combined spectrum is shown later in 
Fig. 3. ft should be noted that the first two 
bins in the SD spectrum were excluded in this 
procedure. We expect these first two bins are 
biased by threshold effects of the order of 10%. 
The deviations of those bins from the Hybrid 
spectrum are in agreement within the system- 
atic uncertainty. 



Figure 1: The energy spectrum multiplied by 
derived from SD using showers at zenith an- 
gles above (filled triangles) and below (opened 
triangles) 60° ([5, 7]), together with the spec- 
trum derived from the hybrid data set (red 
circles) [3]. Arrows indicate 84% CL upper- 
limits [16]. 

Fig. 1 shows the energy spectrum multiplied 
by from SD data using showers at zenith an- 
gles above and below 60° ([5, 7]), together with 
the spectrum derived from the hybrid data set 
(a fluorescence events in coincidence with at 
least one SD station) [3]. The agreement be- 
tween the spectra derived using three different 
methods is good and is underpinned by the 
common method of energy calibration based 
on the FD measurements. Therefore all spec- 
tra are affected by the 22% uncertainty in the 
FD energy scale [1], in which the largest contri- 
bution is the absolute fluorescence yield(14%). 
In this work we have used the fluorescence 
yield reported in [12]. This common uncer- 
tainty does not affect the relative comparison 
of our spectra. The systematic uncertainty in 
the hybrid-only spectrum is dominated by un- 
certainties in the calculation of the exposure 
(16%). The systematic uncertainty in the SD 
spectrum has two contributions: the calcula- 
tion of the exposure (3%) and the statistical 
uncertainty in the calibration of 5'(1000) and 
7Vi9 with the FD energy (<10%). We use 
a maximum likelihood method, together with 
our knowledge of the systematics, to calculate 



The highest end of the spectrum 

Since the 22% systematic uncertainty in the 
energy scale does not modify the shape of the 
spectrum, it is possible to check the continua- 
tion of the spectrum at the highest energies. It 
could be argued that our energy calibration has 
low statistics at the highest energies (see Fig. 
3 in [5]). However no indication of a change in 
the calibration parameters with threshold en- 
ergy used has been found. A dramatic change 
in the hadronic interactions in the energy range 
where no hybrid event is observed could also 
induce false spectral features. However, there 
is no theoretical basis for such a scenario, and 
even if it were the case it will be checked in the 
future with larger statistics in the hybrid data 
set. 

To check the continuation of the spectrum at 
the highest energies we first fit the SD spec- 
trum between lO^^ '^ eV and lO^^-^ eV to a 
power-law function using a binned likelihood 
method. The spectral index obtained is 7 = 
-2.62±0.03(stat)±0.02(sys). The systematic 
error is given by the error on the calibration 
curve in [5]. The number of events expected 
from such a single power-law flux above 10^^'^ 
eV and lO^o eV are 132±9 and 30±2.5 re- 
spectively whereas we observe only 51 events 
and 2 events. Also, the spectral index from 
-[^q19.6 gY i-j-^g highest energy observed 

(1.90±0.16(stat)±0.20(sys))xl020 eV is 7 = 
-4.14±0.42(stat) (Fig.2). A lack of events at 
the highest energies is clear. We then applied 



30th International Cosmic Ray Conference 




2 : — I — I — I — ^ 




18.5 19 19.5 

log(Emin [eV]) 



Figure 2: Upper panel: Spectral index as a 
function of minimum energy in the fit. Lower 
panel: significance (in sigma) of the deviation 
from power-law distribution with spectral in- 
dex from upper panel based on the TP statis- 
tics. 

a statistical test proposed in [13], the so called 
TP-test. The TP statistic allows us to test 
for a power-law distribution on an unbinned 
data set without bias regarding the value of the 
spectral index. Details of this statistical test 
can be found in [14]. The upper panel in Fig. 2 
shows the unbinned maximum-likelihood esti- 
mation of the spectral index (7) and its stan- 
dard deviation (shaded region) as a function 
of minimum energy used in the fitting. A clear 
change of slope at the highest energy can be 
seen. The deviation from the power-law distri- 
bution with 7 shown in this figure is estimated 
based on the TP statistic. The lower panel in 
Fig. 2 shows the estimated deviation in sigma. 
The hypothesis of the pure power-law is then 
rejected with a significance better than 6 sigma 
and 4 sigma for minimum energies of 10^^'^ eV 
and 10^^ eV respectively. 



Astrophysical interpretation 

In the previous scc;tion, we have shown that 
the rejection of the hypothesis of a continua- 
tion of the spectrum in the form of a power-law 
is statistically significant. Moreover, a spectral 
break at ~ 10^^-^ eV, the so-called ankle, is ap- 
parent in Fig. 1. Therefore we fitted the com- 
bined Auger spectrum to the following equa- 
tion: 

J{E-E<Eankle) OC 

J{E-E> Eankle) OC E^^ ^ , ^ . 

l + exp(ii^) 

(1) 

where 71 and 72 are the spectral index before 
and after the break respectively, Eankie is the 
position of the break, and the second term in 
the second equation is a flux suppression term 
where Ec is the energy at which the flux is sup- 
pressed 50% compared to a pure power-law, 
and Wc determines the sharpness of the cut- 
off. Here using a binned likelihood method, the 
values of the parameters obtained are the fol- 
lowing: 71 = -3.30±0.06, 72 = -2.56±0.06, 
Zoffioi;a„fcie-18.65±0.04, logwEc^l9.7A±Qm 
and Wc=0.16±0.04. The xVdof for this fit is 
16.7/16. The black line in Fig. 3 shows the 
result of the fit. 

Fig. 3 shows also a comparison of our data with 

some astrophysical models [9]. These models 
show a flux suppression at the highest energies 
(the GZK steepening [15, 18]). The models ah 
assume an injection spectral index, an expo- 
nential cutoff at an energy of E^ax times the 
charge of the nucleus, and a mass composition 
at the acceleration site as well clS Qb distribution 
of sources. The blue lines in the figure assume 
a mixed composition at the sources, i.e. with 
nuclear abundances similar to those of the low- 
energy galactic cosmic rays. A uniform dis- 
tribution of sources and an injection spectral 
index of -2.2 (close to the shock acceleration 
predictions) are assumed as indicated in the 
figure. Emax is taken as 10^° eV (dashed line) 
and 10^^ eV (solid line). Good agreement is 
found down to energies close to Eankie- Below 
this energy another component is needed. 



Auger Spectrum and Interpretation 



Another set of models which assume only pro- 
ton primaries and Emax = 10^^ eV are shown 
by the red lines. One model assumes uniform 
source distribution with the spectral index - 
2.55 and the other assumes the source evolu- 
tion has a strong redshift dependence (1 + z)^ 
with the spectral index -2.3. It has been sug- 
gested that the spectral break at EanUe can 
be explained as a feature of the propagation of 
a pure proton flux in the extragalactic media 
including pair production [11]. To repro- 
duce our spectrum by this model, we need a 
very stronger source evolution. The distribu- 
tion of the longitudinal profiles of the showers 
observed by the FD also disfavors the pure pro- 
ton assumption [6]. 

Conclusions 

Using data from the southern-hemisphere 
Pierre Auger Observatory, we reject the hy- 
pothesis that the cosmic ray spectrum contin- 
ues in the form of a power-law above an en- 
ergy of 10^^'^ eV with 6 sigma significance. 
This result is independent of the systematic 
uncertainties in the energy scale. A precise 
measurement of the energy spectrum, together 
with anisotropy and mass composition studies 
in this energy range, will shed light on the ori- 
gin of the highest energy particles observed in 
nature. 
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